On the Preparation of Pure States in 
Resonant Microcavities 



Per K. Rekdal^'", Bo-Sture K. Skagerstam^'^''^ and Peter L. Knight^'" 

° QOLS, Blackett Laboratory, Imperial College, London SW7 2BW, 

United Kingdom 

Department of Physics, The Norwegian University of Science and 
Technology, N-7491 Trondheim, Norway 

^ Microtechnology Center at Chalmers MC2, Department of 
Microelectronics and Nanoscience, Chalmers University of Technology 
and Goteborg University, S-412 96, Goteborg, Sweden 

Abstract 

We consider the time evolution of the radiation field (R) and a two-level atom 
(A) in a resonant microcavity in terms of the Jaynes-Cummings model with an 
initial general pure quantum state for the radiation field. It is then shown, us- 
ing the Cauchy-Schwarz inequality and also a Poisson resummation technique, 
that perfect coherence of the atom can in general never be achieved. The atom 
and the radiation field are, however, to a good approximation in a pure state 
\4')A(g)R = \4')a ® in the middle of what has been traditionally called the 
"collapse region" , independent of the initial state of the atoms, provided that 
the initial pure state of the radiation field has a photon number probability 
distribution which is sufficiently peaked and phase differences that do not vary 
significantly around this peak. An approximative analytic expression for the 
quantity Tr[p^(t)], where pA(t) is the reduced density matrix for the atom, 
is derived. We also show that under quite general circumstances an initial 
entangled pure state will be disentangled to the pure state \ip)A(g,R- 
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The Jaynes-Cummings (JC) model [T is one of the simplest but non-trivial examples 
of two interacting quantum systems. It is an important fundamental theoretical model 
of the interaction between a two-level atom and a second-quantized single-mode electro- 
magnetic field (for a review see e.g. Ref.j2j). The model is exactly solvable and realizable 
by experiments e.g. involving the passage of single atoms through a superconducting mi- 
crowave cavity (see e.g. Refs.pj). A direct experimental verification of the one-mode field 
quantization in a microcavity, as assumed in the JC model, has actually been carried out 
Experiments on ion-traps [S] provide for an alternative arena in which case JC-like 
models again are important theoretical tools in describing the relevant physics. Recent 
developments of nano-electronic devices has also led to an experimental realization of cir- 
cuits with superconducting Josephson junctions which behaves as two-level systems 0. 
Entangling such a device to a resonator naturally leads to a description in terms of an 
effective JC model (see e.g. Refs.[7]). The results of the present paper to be discussed 
below may therefore have a quite broad range of potential applications. 

In spite of its apparent simplicity, the JC model has led to many non-trivial and 
unexpected results through the years as e.g. the well-known phenomena of collapses and 
revivals of the atomic population inversion 6|. In Ref. it was found that the atom is 
to a good approximation in a pure state \'iP)a®r in the middle of the "collapse region" 
provided the initial state of the radiation field is a coherent state. It was shown ITU] 
that the existence of this particular pure state can easily be demonstrated in terms of 
the properties of the von Neumann entropy in quantum mechanics. It has also been 
shown that the appearance of this pure atomic state is independent of the form of the 
initial pure atomic state jll . For mixed states of the radiation field it is also possible to 
transfer coherence from atomic states to the state of the radiation field ^21- Here we will 
show that the purification of the atomic state is actually independent of the nature of the 
initial pure state of the radiation field provided that the corresponding number-operator 
probability distribution is sufficiently peaked and that phase differences of this state are 
slowly varying around this peak. We will, however, at this moment restrict ourselves to 
initial states which are not entangled and give some comments on entangled initial states 
in the end of this paper. 

The electromagnetic interaction between a two-level atom, with level separation ujq, 
and a single mode of the radiation field in a cavity with frequency uj is described, in the 
rotating wave approximation, by the JC Hamiltonian jT] (where we take fi= 1) 
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H — uja^a + -uJoGz + g{aa+ + aV_) , (1) 

where the couphng constant g is proportional to the dipole matrix element of the atomic 
transition. Here we make use of the Pauli matrices to describe the two-level atom and 
the notation (t± = {a^ ±icry)/2. The second-quantized single mode electromagnetic field 
is described in a conventional manner by means of an annihilation (creation) operator a 
(a''"), where we have suppressed the cavity mode labels. For g = the atom-field states 
\n, s) = \n) ® \s) are characterized by the quantum numbers n = 0, 1, ... of the oscillator 
and s = ± for the atomic levels with energies En^± = un ±ujq/2. 

For reasons of simplicity we will now consider a resonant system, i.e. uj = ujq. We 
also assume that initially the atom is in the excited state |+) and the radiation field is 
in a general pure state I7) = J2'^=o v^^*""l^)- Hence, the initial state of the system is 
l'^(O)) = 1+) ® I7). The solution to the Schrodinger equation is then 
00 

= ^^e*""e-''*^("+^/2)* cos{gtVnTl)\n,+) -ism{gtVnTl)\n+l,-) ) . 

n=0 

(2) 

In order to describe the evolution of the atom alone it is convenient to introduce the 
reduced density matrix pA{t) — Tr^[ \ilj(t)){ip{t)\ ], where the trace is over a complete set 
of radiation field states. One easily finds that 

PA{t) = +P-{t)\-){-\ + c{t)\+){-\ + cm-){+\ ■ (3) 

Here p+{t) — 1 — p-{t) — Yl'^=o Pn cos^ {gWn + 1) is the well known form for the revival 
probability and 

00 

cit) = -iJ2 VP^iP^cosigtV^^) sin{gtV^)e'^^--^-'^e-'^' . (4) 

n=0 

By straightforward algebra one can now verify that 

Tr[p'^{t)]=Tr[pl{t)]=pl{t)+pUt) + 2\c{t)\' , (5) 

where the radiation field density matrix is P'y{t) — Tr^[ \ip{t)) {ip{t) \ ], and where the trace 
is over the atomic states. This means that purity occurs in both subsystems at the same 
time and precisely at the same rate. We notice that, in general, the cu dependence cancels 
out in Tr[ p^(t) ]. For a coherent state 1^;), with = exp(iQ;), we now have that = na 
with photon-numbers that are Poisson distributed, i.e. p„ = e-xp{—\z\'^)\z\'^'^ /n\. The a 
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dependence will therefore also cancel out in Tr[ p^(t) ] for any initial state I7) which 
has the same phase dependence as a coherent state. By a tedious but straightforward 
calculation, the identity Tr[ p^{t) ] = Tr[ p\{t) ] can actually be shown to be valid for any 
initial pure state of the form [^^(O)) = + b\—)) ® I7) with |ap + |6p = 1. 
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Figure 1: This figure shows T\:[p\{t) ] (= Tr[p^(t) ]) as a function of gt with an initial coherent 
state of the radiation field with |ap = 49 based on the exact numerical results due to Eq. 
This equation with a Gaussian photon- number distribution p„, with mean n = |ap and standard 
deviation a-a = \/^) describes the data with a very high numerical accuracy. With our choice 
of parameters and within the accuracy of the figure shown, Eq. © and the approximative 
expression of Tr[p^(t)] according to Eqs. p4|) and agree. As the mean number of photons 
gets larger Tr[p^(t ~ to = ^V^)] approaches one at an exponential rate as can be seen from 
Eq. (UHl). 

The purity of the atomic state can then be determined considering the time evolution 
of the quantity Tr[ p\{t) ][9 . A necessary and sufficient condition for the ensemble 
to be described in terms of a pure state is then that Tr[ p^(t) ] = 1, in which case 
clearly a state-vector description of each individual system of the ensemble is possible. 
On the other hand, for a two-level system, a maximally mixed ensemble corresponds to 
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Tr[p\{t) ] = 1/2. Alternatively one may consider the von Neumann entropy of the density 
matrix pA{t) 10 or simply its eigenvalues JI]. One can also make use of the Schmidt 
decomposition of composite system (see e.g. Ref.[T!^) to reach the same results as given 
below. 

The definition of c{t) according to Eq. (0)) above now suggests an immediate interpreta- 
tion in terms of an I2 Hilbert space scalar product of two complex vectors a and b with com- 
ponents an = ^Jp^co?,{gt\/n + 1) exp(ia„) and 6„ = ^"Pn-x sm{gty/n) exp(ia„_i) respec- 
tively with bo = 0. The Cauchy-Schwarz inequality then tells us that |c(t)p < p+(t)p„(t) 
with equality if and only if the vectors a and b are parallel, i.e. a„ = j3{t)bn where 
= p+(t)/j9_(t). In general we then see that 1/2 < Tr[ p^(t) ] < 1. The condition 
a„ = P(t)bn can be used to find a particular pn such that Tr[ p\(t) ] = 1 at a fixed time 
t = tf. In order to see this let us assume that a„ = an + a^. At t = tf we then obtain 
a recursion formula p„ cos^i^gtfy/n + 1) = \P{tf)\'^pn-ism^{gtfy/n). Depending on the 
parameter gtf this recursion formula can be solved. The corresponding probability distri- 
bution can now be used as an initial distribution and is then such that Tt[ p\{t = tf) ] = 1 
by construction. The time parameter tf is, however, not in general related to any natural 
revival time of the system. 

In order to proceed in a more general setting, we assume that the distribution pn, 
with mean n and variance (T„, is peaked around n = n or more precisely cr„/fi ^ 1. We 
observe that cr„ should, however, not be arbitrarily small in order to have a non-zero 
c{t). An explicit condition on (j„ will be discussed below using Poisson resummation 
techniques. Following the analysis of Ref. P, and by considering time-scales t close to 
half of the first revival time trev = lT^\Jn + 1/(7, i.e. t 'r^ t^ = trev/'^i one finds that 
Tr[p^(t) ] = l/2 + 2|c(t)p using that p±{t ~ to) = 1/2. As long as ^|t-to| << ^^/fiTl, 
we can make use of the fact that gt\/n + 1 = gt^fn^i^ jl [Oj and we find that |c(t)p = 1/4 
provided Y^exp(ia„) = Py/Pn-i exp(m„_i) for all n with a complex phase P = 1). 
This is possible with high accuracy only if i) pn is non-zero (and constant) for a finite range 
of n as we actually have assumed, and ii) that the phases a„ are of the form a„ = an + ao, 
i.e. the phase differences a„ — a;„_i do not vary significantly around n = n. Under such 
circumstances we then see that Tr[ p\(t) ] = 1 and the atom is in a pure state at t ~ to- 

In FigH we illustrate the behaviour of Tr[p^(t)] in the case of an initial coherent 
state for the radiation field with n ^ 1. The distribution pn can the be approximated 
with a sharp Gaussian distribution obeying the condition i) above with high accuracy. 
In Fig. 121 we illustrate the behaviour of Tr[p^(t)] in the case of a pure state with pn 
which varies rapidly in terms of a Schrodinger cat state here chosen to be of the form 
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|7) = N{\a) +exp{i(j))\ - a)) with \N\^ = 1/(2 + 2 cos(0) exp(-2|an (see e.g. Ref.|ISI). 



For such a state it follows from the analysis above that \c{t)\ cannot reach its maximal 
value 1/2 as is also clearly exhibited in Fig. |21 The behaviour of Schrodinger cat states 
at t = to including the effects of damping has been studied in more detail in Ref. |16j. 

It is now clear from the definition Eq. Q of c{t) that if the phases an are rapidly 
varying functions of n, c{t) will be small or zero. If in particular the phases are random 
the state of the radiation field will correspond to a mixture instead of a pure state and 
we have to average over the phases of pA{t) before computing its square. We then obtain 
Tr[p^(i(:)] = 1 + 2p^(t){p^(t) — 1) with purity of the atomic state only at times t such 
that p+{t) = 1 or 0. In such a case we also notice that Tr[p^(t)] in general is different 
from Tr[p^(t) ]. The purity of the state of the radiation field is therefore required in order 
to purify a general initial mixture of the atom. For related studies of purification in the 
JC model see Refs.[T7j. 

The quantity Tr[p^(t)] can now be re-written in a form where quantum revivals 
become explicitly by making use of a Poisson summation technique ITl!j. As above, we 
consider a probability distribution p„ is peaked around n = n with cr„/n <^ 1. If the 
phase differences a„ — are slowly varying functions of n around n = ^ 1, i.e. we 
assume that ~ an + ao for n n, we then see that \c{t)\ according to Eq. (jl} can be 
written in the form 




(6) 



n=0 



and we can then write Tr [ p\ (t) ] in the following convenient form 



(7) 



where we have defined the functions 



oo 



oo 
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n=0 



and 



oo 




(9) 



n=0 

An exact Poisson resummation technique ^1] now enables us to e.g. write 



oo 



,2igt 




(10) 
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Figure 2: The quantity Tr[p^(t)] (= Tr[p^(f)]) as a function of gt when the initial state 
of the radiation field is a Schrodinger cat state where pn = lap" [1 + (— l)"cos</>] /[n!(el"l^ + 
e-l"P cos(/)) ] with |ap = 49 and = 0. Purity of the atom occurs in this case for gt = only. 

where 
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Here p{n) is an analytical continuation of p„ and Sp{n) = nun — gt\Jn + 1. If p{n) 
varies slowly as compared to the variation of Sy{n) one can evaluate f^it) for z/ 7^ by 
a stationary phase approximation [11]. For a Gaussian distribution with mean n and 
variance a^, i.e. p{n) = exp(— — 72)^/20-^)/ y^27raf, this corresponds to the condition 
cr^ ^ (n + l)/7rz/, provided n ^ 1 . The contribution from v = Q can be estimated by 
again considering a Gaussian distribution with n ^ 1 in which case one finds that 



dn p{n)e 



~ Q-igto-nY /2n^2igtVn+l 



(12) 



The stationary phase condition on Sv{n) now immediately leads to the general revival 
times trev = ^Tcuy/n + l/g. Applying the same techniques to S(t) according to Eq. Q 
leads to revival times which are of the order n larger than the those of Sc{t) and Ss{t) 
and we therefore make the approximation 
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Figure 3: This figure sliows illustrates the effect of squeezing on Tr[/9^(t)] (= Tr[p^(t)]) as 
a function of gt when the initial state of the radiation field is a squeezed coherent state I7) = 
S{r)\a) with r = 0.75 and a real a (~ 14.72) chosen such that the mean number of photons 
n = 49 is the same as in Fig. ^ The photon-number distribution is well approximated by a 
"squeezed" Gaussian distribution with mean n and o"^ = n^'^^. 



n=0 



sm 



gt 



dn p{n) sin 



gt 



g -(gto„)V32n3 ^.^ / gt 



(13) 

where we again have made use of Gaussian distribution and n ^ 1. The approximations 
Eqs.(jni), and (fT^ together with the stationary approximation above enables us to 
write the quantity Tr[ p\{t) ] in an analytical form, i.e. 

gt 



Tr[ p\{t) 



1 1 

2 + 



where 



Wo{t) 



and for u > 1 



g -(gto„)2/fi _ 2 g-(gt(T„)2(l+l/16n2)/2n / gt 



sm {2gtVn + 1) 



Wu{t) 



gt 



IT ■ 



V 



■ gt 

■2nu 



o 9t ( f gt ^2 

— 2 ; p' ' ^ 



= sm 

1/ V27rt/ 4 



(14) 



(15) 



(16) 



Apart from the approximations used in order to obtain Eqs. ()12j) - (fT3j) . the Poisson re- 
summed expression for Tr[ p\(t) ] as given by Eq. with Eqs. (fTS|) - (fT^ is vahd for any 
distribution p{n) provided that n ^ 1 and 

(n+ l)/z/7r < (T^ < . (17) 

With this condition, which is always satisfied for a Poissonian distribution, we have 
also verified through numerical calculations that the various approximations leading to 
Eq. (jl4|) describe a much larger class of p„-distributions than just Gaussian ones. For 
n ^ 1 we can now extract the leading behaviour of Tr[ p\{t) ] close to t = trevf^, i-e. 
11 

which show the exponential approach to purity. These results shows that the quantity 
Tr[ p\{t) ] exhibits revivals not only for an initial coherent state, as one is lead to believe 
from Ref. but for any initial state of the radiation field in the form of a pure state 
I7) — Yl'^=o Vp^^^°'"\^) "with a sufficiently peaked probability distribution p„, satisfying 
Eq. (fT7|). and with phase differences a„ — varying slowly around this peak. 

In particular, if the radiation field is initially in a coherent state with average photon 
number n = |ap as in Ref. [HI, then p„ is Poisson distributed and reduces to the Gaussian 
distribution with variance = n if n ^ 1. In this case, Eq. p6|) can be reduced to 
1 r {gty (gt-gty)^ 



Tr[ plit ^ hre.)] + e-'^"'^'/^^^' + e-'^"-' - 26-'^"-'/^ sin(27rn)) , (18) 



e 



- ^ e e — sin ( ^ ) sin ( ^ ) , 19 

where 

gU = 2nuy/fi . (20) 

As shown in Fig. Qthe analytical form Eq. (|14p with Eqs. (jl5|) and (|19|) describes the actual 
from of Tr[ p\{t) ] with a very high numerical accuracy. If we make the photon- number 
distribution p„ more peaked as compared to e.g. a Poissonian distribution we now see 
from Eq. (fT^ that the approach to purity at later times tren/S + 27r/c, with k = 1,2, 
will be more visible. In Fig. |21 we illustrate this feature for a squeezed coherent state 
I7) = S{r)\a) with the same mean-value of photons as in Fig. ^ Here we consider a 
squeezing operator S{r) = exp(r(a^ — a"'"^)/2) and a real leading to the photon-number 
distribution [TS] 

tanh(r)'"' / 2/ , / \\\ tj f ^ ^ ^ 



Vn 



, V, exp f - |ap(l -tanh(r))WJ^^^) , (21) 
n!2"cosh(r) ' ^ "V^sinh(2r)/ ^ ^ 
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where Hn is a Hermite polynomial. In this case the Poisson resummed expression Eq. ()14|) 
reproduces the exact answer, within the numerical accuracy of Fig. El by making use 
of the fact that the photon- number distribution is well approximated by a "squeezed" 
Gaussian distribution with mean n and o"^ = ttP'^^ except for times close to to- For the 
parameters chosen the inequality Eq. (fT7j) actually fails ioi v = 1. Evaluating the integral 
Eq. more carefully for u = 1, which e.g. can be done analytically by making use of 
a Gaussian approximation for p„, one restores the agreement between the exact answer 
and the approximation Eq. (fT^ with high numerical accuracy. 

In concluding we would like to comment on the nature of pure state \iP)a(^r at t ^ 
^0 = trevf^ = irVn + !/(?. We again assume that n ^ 1 and that the conditions Eq. (fTTj) 
are fulfilled for any probability distribution pn under consideration. As an initial state we 
consider a general pure entangled state, i.e. 

|V;)=a|+)®|7+) + 6|-)®|7-) , (22) 

where |ap + |6p = 1 and |7±) = '\/p^<2*"" \n). As long as \t — to\ « l^ft^T^lg we 

again make use of the fact that gt^n + 1 = gt-Jn + 7r/2 P and we find that 

|^(t ^ to)) ^ (1+) + |-)e^-*°e-^^") ® |7(t ^ to)) , (23) 
with Aa = Aq;^(?7,), where IS.oc^{n) = — a^_i , and 

oo 

l7('^)) = - E e-^-("+^/')* (ae^"" v^sin(^tv^) + ^6e^-"+i ^p~^, co^{gt^) ) \n) . (24) 

n=0 

In order to obtain Eq. (j23|) we only assumed that the distributions have their main 
contribution at the same n = n around which Aa^{n) are assumed to be slowly varying 
and equal when evaluated at n = n. Apart from the phase factor e*'^*'', Eq. ()23|) agrees with 
the result of Ref.0. With regard to experimental realizations it is, however, important 
to reahze that the relative phase in Eq. depends on both Aa and uto. The JC- 
model therefore predicts disentanglement at t ~ trevf^ to the pure atomic state \4')a = 
(1+) + |— )e*'^*''e~*^")/-\/2 fairly independent of the nature of the initial pure entangled 
state. 

In summary, we have shown that exact coherence of the atom is in general never 
regained for a JC model with a general initial pure quantum state I7) of the radiation 
field. One can, however, come arbitrarily close to a disentangled pure atomic state \iP)a(^r 
when I7) has a general form ^^0 y/PnG^°'"\n) provided that the probability distribution 
Pn is sufficiently peaked around its mean value n ^ 1 and that the phases differences 
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«n — do not vary significantly around this peak. Under such conditions we have 
then derived an approximative analytical expression for the quantity Tr[ p\{t) ]. Hence, 
the quantity Tr[ p\{t) ] exhibits revivals not only for an initial coherent state of the 
radiation field, as one is lead to believe by the results of Ref . j2| , but actually for a very 
general set of pures of the radiation field. By making use of the analysis of Ref. [TJI one 
can, in a straightforward manner, make use of the methods above and verify that this 
actually remains true independent of the initial state of the atom. We have also seen that 
under quite general circumstances an initial pure entangled state also leads to the same 
purification of the atomic state. 
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